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O" 1 Abstract 

In this paper, we investigate the nature of self-interacting potential 
that can support the accelerated expansion of the universe as indicated 
by observational data. For this purpose, we consider the Bianchi type 
I (BI) universe model in the Brans-Dicke (BD) field scenario and use 
a power law for both the scalar field and the scale factor. It is found 
that the accelerated expansion of the universe can be discussed for a 
positive power law potential with negative values of the BD parameter. 
We also explore the evolution of energy density perturbations. This 
analysis indicates that the model allows growing modes for negative 
values of the BD parameter and m > 0. 
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1 Introduction 

Many astronomical experiments such as Supernova (la) 1,2 \ Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) 3 \ Sloan Digital Sky Survey (SDSS) 4 \ 
galactic cluster emission of X-rays % large scale-structure 6 ) and weak lensing 
7 ) , provide evidence of an accelerated expansion of the universe in the present 



* Email address: msharif.math@pu.edu.pk 
tsmathematics@hotmail.coni 



1 



phase. As the baryonic matter and dark matter do not provide the total en- 
ergy density of the universe, it has been suggested that there is some kind 
of " missing energy" . The Chaplygin gas, phantom, quintessence and cosmo- 
logical constant 8 ' 9 ^ are some appropriate candidates of dark energy (missing 
energy), considered to be responsible for the expansion of the universe. In 
order to conciliate the nucleosynthesis scenario in the radiation-dominated 
era, it is argued that the universe was expanding with deceleration in its 
early phase. 

Quintessence is a dynamical slowly evolving spatially inhomogeneous com- 
ponent of energy density with a negative pressure 9 ) having an EoS parameter 
lying between and —1. However, this model was proved to be not very suc- 
cessful as it faces problems like fine tuning and cosmic coincidence 10 ). To 
solve these problems, a special form of quintessence field called the tracker 
field (with a wide range of initial values of p^, such fields roll down along 
a common evolutionary track with the present energy density of the uni- 
verse p n and end up in the observable universe with p^ comparable to p n 
in the present epoch) has been constructed u \ Different types of potential 
have been developed 12 ) to understand the nature of dark energy using a 
quintessence model but each of them has its own drawbacks. 

Scalar-tensor theories of gravity as modified gravity theories are proved 
to be very effective in explaining many outstanding problems in cosmology 
13-16) jjkg i n f| a tion, the "early- and late-" time behaviors of the universe, 
coincidence problem and cosmic acceleration. Brans and Dicke proposed a 
scalar-tensor generalization of general relativity (GR) by introducing a time- 
varying gravitational constant (as G(t) = ^j) and the direct interaction 
of a scalar field with a geometry known as the BD theory of gravity. The 
weak equivalence principle, Mach's principle 1S ) and Dirac's large number 
hypothesis are also accommodated in this theory. Using local gravity 
tests, it has been shown 20 ' 21 ) that the generic dimensionless BD parameter 
uj should be very large, i.e., u > 40,000, for being compatible with solar 
system experiment bounds. This theory goes over to GR if uu — > oo and the 
scalar field is constant 22 \ In general, it has been shown 23 ' 24 ) that the BD 
theory is reducible to GR only if the trace of the energy-momentum tensor 
does not vanish, i.e., ^ 0. 

Many researchers have explored the cosmic evolution of the universe in 
the BD theory. In the GR and scalar tensor theories, much work has been 
done using Bianchi models to discuss different stages of the cosmic evolution 
25,26) _ j n re f 27 ; it has been pointed out that some large-angle anomalies are 
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seen in CMB radiations that violate the statistical isotropy of the observable 
universe. For a better description of these anomalies, plane symmetric and 
homogeneous but anisotropic universe models play a very significant role. 
Moreover, it is pointed out 28 - ) that removing a Bianchi component from 
WMAP data can explain various large-angle anomalies yielding an isotropic 
universe. Thus, the universe may have accomplished a slight anisotropic 
geometry in cosmological models irrespective of inflation. King and Coles 
29 -* described the evolution of magnetized axisymmetric Bianchi I universe 
in the presence of vacuum energy. Singh and Rai 30 ) constructed various 
BD cosmological models using Bianchi models with a perfect fluid to explain 
the evolution of the universe. In our recent study 31 \ we have investigated 
some Bianchi type I universe models using perfect, anisotropic and magne- 
tized anisotropic fluids in self-interacting BD cosmology. Bermann 32 ^ investi- 
gated different cosmological models using a constant deceleration parameter. 
Khoury and Weltman 33 ) explored the coupling of a scalar field with matter 
of unity order, which is allowed to evolve cosmologically. 

Bertolami and Martins 13 ^ obtained accelerated expanding solutions using 
quadratic potentials. They also evaluated energy density perturbations and 
discussed the growing modes of their model. By extending the concept of 
quintessence to flat non minimally coupled scalar tensor theories, Ritis et 
al. 34 ) found a family of exact solutions parameterized by the parameter V. 
These solutions are proposed as a class of tracker fields known as tracking 
solutions. Sen and Sen 35 ) investigated accelerating solutions using a potential 
with a time- dependent mass squared term in BD cosmology. By introducing 
the non minimal coupling term ipR^-, different kinds of potential have been 
discussed 36 ) . Sen and Seshadri 37 ) investigated the nature of self- interacting 
potential as well as the evolution of energy density perturbations in this 
theory. 

In this study, we explore the nature of self-interacting potential and BD 
field energy density as dark energy, which can account for the accelerated 
expansion of the universe by taking a perfect fluid and the LRS Bianchi type 
I (BI) universe into account. This paper is organized as follows. In the next 
section, we formulate the field equations of the self-interacting BD theory for 
the BI universe. We discuss solutions by taking a power law ansatz. In §3, we 
provide the evolution of energy density perturbations by taking a perturbed 
form of field equations. We discuss the results in the last section. 
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2 Self-Interacting Potential and Dark Energy 

The action for the BD theory with the self-interacting potential V(4>) is given 

by 37) 

S = J d^x^g^R - j<P ,a <t>,a ~ V{4>) + L m ], a = 0, 1, 2, 3. (1) 

Here, L m is the matter part of the Lagrangian and uq is the constant BD 
coupling parameter (where the units 8ttGo = c = 1 are imposed). The 
corresponding field equations derived from the action (pQ) are 



U ° \A A 1 A AM j_ 1 \A rul J_ T t" J_ KM 

- 2^0,«0 J + - ^°0J + — + -g^-t 



The terms T = g^T^, □ = A M A M and A M represent the trace of the energy- 
momentum tensor, the d'Alembertian operator, and the covariant derivative 
respectively. 

The homogeneous, anisotropic, and spatially flat BI universe model is 
given by 38 - ) 

ds 2 = dt 2 -A 2 (t)dx 2 -B 2 (t)(dy 2 + dz 2 ), (4) 

where A and B are the scale factors along the transverse direction x and 
the two equivalent longitudinal directions y and z, respectively. The energy- 
momentum tensor for a perfect fluid is 

Tfo, = {p + P)u^u„ - Pg^, (5) 

where p and P are the energy density and pressure, respectively, and w M is 
the four- velocity with the normalization condition u^u 11 = 1. The energy 
conservation equation is 

P+(^ + 2§)(p + P) = 0. (6) 
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The corresponding BD field equations become 



AB B 2 2 2 A B } (f) 20 ' 1 ' 



B B 2 2 2 5 20 ' 1 ' 

B , i ^ P q; o 2 A B j> V{4>) 

B A AB 2 2 M 5 j 20' U 

while the Klein-Gordan equation for scalar field is 

The average scale factor and mean Hubble parameter respectively turn out 
to be 

c?(t) = AB\ H(t) = ±(± + 2 B I ), 
whereas the expansion and shear scalar are 

a A ~B I .A B. 

It is reported in ref. 39 that, for a spatially homogeneous universe model, 
the normal congruence to homogeneous expansion yields the ratio | as a 
constant, i.e., "the expansion scalar 9 is proportional to the shear scalar 
a" . This physical condition leads to the following relation between the scale 
factors: 

A = B m , 

where m ^ 1 is any constant (for m — 1, it reduces to a flat FRW model). 
This condition has been used to investigate various exact universe models 
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40 42) ^his assumption leads to the following set of BD dynamical equations: 

9^ + ^! _ 2 H_l + YM ( 12 ) 

B 5 2 2 2 B 20 ' 1 J 

(m + ij 5 +m J B 2 " 2 2 (m+ij 50 + 20 ' [L6) 



p = _( m + 2 )|(p + P). (15) 

Here, the unknown functions are B, V, p, P and 0, and there are only three 
independent field equations. In order to have a closed system of equations, 
we take the power laws for the scale factor and scalar field as 

B = b (lr, = o (f) /3 , (16) 

where t , b , and O indicate the present time, scale factor and scalar field, 
respectively. The deceleration parameter is given by 

q=-(l 



aim + 2) 



For an accelerated expanding solution of the field equations, the param- 
eter a must be greater than ID ^ L - Adding eqs. (fl2~]) and ( |T3l) . it follows that 



(m + 3)- + (m +l)_ = - Wo (-) -( m + 3)---2--2- + —j—. (17) 



Equation (|T5l) leads to 

P + P = (18) 

{m + 2)B K J 

Subtracting eq. ljTTI) from eq. (fT71) and substituting the power law ansatz given 
by eq. (|T6|) along with eq.( !T8|) in the resulting equation, the energy density 
for matter can be written as 

P = Pc {m) t p ~\ (19) 
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where 

(m + 2)(p a (m 2 - 3m + 2)a 2 /3(ma + a + 2) (m + 3)a 



(/3-2)4 1 2 2 2 

+ (l + co )(3 2 }. (20) 



Pc{m) 



Using the energy density given by eq. ( 1191) in eq. ( 1151) , the corresponding pres- 
sure becomes 

P = Pcim)^- 2 , (21) 

where 

0o(2 -p-(m + 2)a) Am 2 - 3m + 2)a 2 _ f3{ma + a + 2) 
c(m) ~ {^2)4 1 2 2 



(m + 3)a 2 
+ (1 +u )(J \. 
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From these expressions, we can write EoS as 

P = 1BP, 1b= f~f , -l, (22) 
a (to + 2) 

where jb is the EoS parameter. Since the fluid under consideration is a 
perfect fluid (0 < 7b < 1), it constrains the parameter (3 as 

2 - 2a(m + 2) < 8 < 2 - a(m + 2). (23) 

We would like to mention here that for any positive m other than 1 and 
a > — 777, 8 remains negative. 

m+2 ' " o 

Solving the field equations ( TTT]) - (|T4|) for the self-interacting potential, we 
obtain 

C/3-2) 

V{4>) = V c <j)—, (24) 



where 



i /|3 r , r m 2 + 5m+6 ,2a(m + 2) - 2 + 8 w m 2 - 3m + 2., 

v < = ir [a< — 2 < — )< — s — » 



o r^2 



/ m + 3,(2fl + 2a(m + 2)-4) „ r 3m + 7 a m+la+2 

+ a( )— — - — - + 1 + - 

K 2 ' 2-/3 ,x 2 2 

2a(m + 2)-2 + /3 1 2 (1 + u; )(2^(m + 2) - 2 + g) 

X (2 = 0) i + (2 = 0) ^ 



Now, we discuss the energy density and pressure for the missing energy (dark 
energy) by taking the energy density and pressure due to the scalar field. 
These are given by 

Uof_ V (m + 2)jj 
H ~ YJ + 1 2 B ' 

P * = YJ'2 +<P+ 2B ■ 
After inserting their respective values, these quantities become 

P<t> = Pc W t?- 2 , P<p = P cW t fi ~ 2 , (26) 

where 

PcW = [-J a/3(m + 2) + -(a { ( ) 

2a(m + 2) - 2 + f3 (m + 3) (2/3 + 2a(m + 2) - 4) 

X ^ 2-/3 ^ + a 2 2-/3 

(3m + 7)a ((m + l)a + 2) (2a(m + 2) - 2 + g) 

+ l{ 2 + 2 (2-/3) 1 

R 2 h _ (l + ^ )(2q(m + 2)-2 + /3) O (t) -, 
+ ' { (2-/3) jjJ ^' [ V 



u l3 2 (m + 3) 2 1, 2r ^ 2 + 5m + 6 m 2 - 3m 
^eW = [— + a/3 v 2 + /3 2 -/3--(a 2 { ( 

lW 2a(m + 2) -2 + /3 xl (m + 3) (2/3 + 2a(m + 2) - 4) 

+ 1)( 2—]3 )} + a — 2^ 

(3m + 7)a ((m + l)q + 2) (2a(m + 2) - 2 + /3) 

^ { 2 2 (2-/3) 1 



.2/, (l + ^ )(2q(m + 2)-2 + /3) O 



The energy density and pressure for the scalar field obey the EoS, i.e., P^ = 
70P^, where 

= _ t + w + Mo) - g _^ a g (m) 
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In order to have a positive energy density for both matter and scalar field, 
the BD parameter should be constrained as 

2(2-/3) . „. . (m 2 + 5m + 6)a 2 ,m 2 -3m 



/3 2 (4-2/3-2a(m + 2)) L ^ v ' 4 v 4 

2 a 2 (2a(m + 2) +0-2) a(m + 3)(2a(m + 2) + 2/3 - 4) 
Xa + 1P 2^ 4(2-/3) 

/3(3m + 7)a /3 /3(2a(m + 2) + /3 - 2)(2 + a(m + 1)) /3 2 



4 2 4(2 -/3) 2 

/3 2 (2a(m + 2) -2 + /3) 2 + a(m + l) a(m + 3) 

< ^0 < [ tt; 1 1 



2(2-/3) J 2/3 2/3 2 

a 2 (m 2 - 3m + 2) 



2/?2 (30) 

The deceleration parameter requires a > for an accelerated expansion 
of the universe. Since m > 0, (m / 1), we can discuss two possible values of 
m, i.e., m > 1 and < m < 1. The feasible region for uo allowed by eq. fl30|) 
as (/3, oj q ) space is shown in Fig. 1. Fig. 1(a) indicates the (/3, u ) space 
for m > 1 and Fig. 1(b) shows the allowed region for these parameters 
with < m < 1. These figures indicate that the positivity condition for 
both energy densities implies negative values of the BD parameter in both 
cases. We found that the admissible range of BD parameter turns out to 
be —1.5 < uq < — 1 only for < m < 1 with /3 < —5.5. This range for 
the BD parameter is negative. However, it does not yield ghost instabilities. 
A ghost cannot exist consistently, however, in the present model, the ghost 
field could be avoided by a suitable selection of parameters. 

Using eq.fllip. the matter and scalar field density parameters are defined 
as 2 2 

Vt m = — : — , VLa, = — : — . (31) 

<f){l + 2m)B 2 <f){l + 2m)B 2 V ; 

Inserting B and from eq. (I16|) . it follows that 

;i + 2m)a 2 ' (l + 2m)a 2 

In terms of matter and the scalar density parameters fl m and fl^, the EoS 
parameter 7^ in eq.( j25]) can be rewritten as 

q/3(l+m) _ (P-2)tf ) a(l+2m)n m , (m 2 -3m+2)tt 2 _ (m+3)a 

^ = ^ (1 + 2m)a 2 % • (33) 
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(a) feasible region for o (b) feasible re S ion for <*> 

(i)q axis 




/3 axis 



Figure 1: Admissible region for u vs (3 for (a) m = 2 and a = 0.8, (b) 
m = 3/4 and a = 1.1 (Color online). 



(a) feasible {a,j}) space ( b ) feasible space 




Figure 2: Allowed region for parameter (3 versus a with (a) 7^ = —1.24, m = 
2 and (b) 7^ = —0.6, m = 3/4 (Color online). 



The observational data provided by the Supernova Cosmology Project 
and High z Survey Project can lead to an accelerated expanding universe 
with a positive cosmological constant By taking the best-fit values of 
the density parameters for matter and the scalar field (for flatness of the 
universe, i.e., Vt m = 0.28, Vl^ = 0.72 ) 37 ^ as well as the bounds for (5 given by 
eq. (!23|) . we can determine the range for the EoS parameter 7^. We consider 
the case when < m < 1 and take m = 3/4. The corresponding constraints 
for the other parameters are a > 1.09 and —8.45 < (3 < —3.225. Using these 
values, the constraint for the EoS parameter 7^ turns out to be —1.44 < 
70 < 0.71. Likewise, in the case of m > 1, we take m = 2 and thus we 
have a > 0.75 and —4.4 < < —1.2. The corresponding bounds for 7^ 
will be —1.8 < 70 < —0.23. The allowed region in the (a, 0) space with 
0.2 < Q m < 0.4, which satisfies -1.44 < 7^ < 0.71 and -1.8 < 70 < -0.23, 
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-4 



/3 axis 



Figure 3: (a) Plot of H t vs (3 for a = 0.8, Q m = 0.28, and m = 2. Here, 
the pink horizontal line corresponds to a = 0.8. (b) Plot of Hoto vs /3 for 
a = 1.1, Q m = 0.28, and m = 3/4. Here, the pink horizontal line indicates 
a — 1.1 (Color online). 

is shown as a shaded region in Fig. 2. In Fig. 2(a), we take the case m > 1, 
and Fig. 2(b) corresponds to the case < m < 1. 

As the BD theory of gravity has a salient feature of a time-varying grav- 
itational constant, i.e., G(t) = ^y, the rate of change in the gravitational 
constant at the present time can be calculated as 



where H is the directional Hubble parameter at the present time and we 
have used the power law ansatz eq.( !T6|) . For a viable model representing the 
accelerated expansion of the universe, the constraint [§]t=t < 10~ 10 should 
be satisfied 39 ^ . For our model, the best-fit value of the parameter a = H t 
admitting the accelerated expansion can be written as 



The graphical behavior of this expression in the cases < m < 1 and m > 1 
along with appropriate values of other parameters is shown in Fig. 3. It can 
be seen that the graphs of H t with the allowed range of u intersect at least 
at one point on the respective horizontal lines a = 0.8 and a = 1.1 in both 
cases. 

The correction to Newton's law leads to two kind of problems: the varia- 
tion in the gravitational constant and the propagation of the fifth force. From 




(34) 



(#oto) 2 



a(m + 2) m 2 — 3m + 2 2 (3(ma + a + 2) (m + 3)a 



(l + 2m)Q m l 2 2 2 

+ {l + ^n (35) 
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the cosmological point of view, the gravitational constant should be positive 
(it restricts 0o to be positive) and its variation should be constrained within 
certain limits 13 ~ 15 ). The case of a constant BD parameter yields small vari- 
ations that should be constrained by [^]o < 4 x 10~ 10 yrs. In our case, the 

expression is given by [^] = — where the present time is taken to be 
t = 14 ± 2 Gyr. By choosing the appropriate values of the parameters, it 
can be shown that in our case this constraint can be easily satisfied. For 
example, if we take a = 0.8, m = 2 and /3 = —1.3 as —4.4 < (3 < —1.2, then 
we have [^] = 0.93 x 10~ 10 , which safely lies within the limit < 4 x 10~ 10 . 
Thus it can be concluded that, in our case, the variation in the gravitational 
constant (due to presence of scalar field) could satisfy the observational limit 
suggested for cosmic acceleration. 

The Brans-Dicke theory in the Jordan frame can be transformed to the 
Einstein frame by defining a conformal transformation. The resulting action 
in the Einstein frame corresponds to a coupled quintessence field scenario 
with the coupling given by 2Q 2 = 2aJB ^ +3 (this leads to a constant coupling 
Q = — in metric f(R) gravity) 43 ~ 47 ). It is argued that when this cou- 
pling is of order unity, local gravity tests become inconsistent owing to the 
propagation of fifth force between the field and the non relativistic matter. 
Whenever the gravity is modified from GR, there are some constraints that 
come from the local gravity tests like solar system tests, and the violation of 
the equivalence principle etc. 

The action of the BD theory in the Jordan frame can be transformed 
to the Einstein frame via conformal transformation defined by the coupling 
factor F = e~ 2Q< ^ (basically, we define g^ u = F{<j})g llv ). For a massless scalar 
field, the coupling factor should be constrained by the condition \Q\ < 2.5 x 
10~ 3 that comes from the experimental bound on the BD parameter uj BD > 
40, 000. When the field is massive, this coupling factor can be large \Q\ ~ 1, 
consequently, the Chameleon mechanism provides a way to be consistent with 
local gravity tests. In this regard, the field mass needs to be small in order to 
realize the present acceleration on cosmological scales pa H , while it is 
large in the region of high density to avoid the propagation of the fifth force 
^> H . The field that changes its mass depending on the environment 
is the Chameleon field. In a spherically symmetric configuration, the field 
equation in the Einstein frame is given by 

d 2 <P | 2d<P dUeff _ Q 

dr 2 r dr d(p 
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where U e ff — U + e Q< ^p m . A detailed discussion on thick- and thin-shell 
solutions in this context by finding the minima of the effective potential is 
available in the literature 43_47 ). 

The thin-shell parameter is defined to be e t h = ^qq^ a an d its upper bound 
can be found using solar system tests or the violation of the equivalence prin- 
ciple. If the thin-shell parameter is much less than 1, i.e., e th < 1,9 q* — , then 
solar system tests are satisfied even for \Q\ = 0(1). Likewise, the propa- 
gation of the fifth force (defined by F = —6Q 2 £th^y^) will be suppressed 
as compared to gravitational force if eth < 8-8 q2° ? , which further leads to 
the constraint \<ps\ < 3.7 x 10 -15 , where \4>b\ \4>a\- This analysis is 
then applied to a runaway form of the potential, e.g., the inverse power 
law potential, and the respective constraint on M is discussed 45 > 46 ). In our 
case, we have found the self interacting potential with a positive power as 

,8-2 a_2 

V = V c (f) (* , > 0, for negative p values. In the present case, the minima 
of the effective potential inside and outside a spherically symmetric body 
(Pa > Pb) are given by 



<Pa = <Pb = 

QpA QpB 

In order to suppress the propagation of the fifth force, we must have — nV c < 
(3.7x 10 _15 ) 1_n <5p_B, where ps can be taken as 10~ 24 g/cm 3 . Moreover, for the 
Chameleon mechanism to work, the mass of the field should be constrained 
by the condition itla > Hq (this allows the condition 1/(itlaTc) «C 1) which 

2-n 

in our case results — j^fyj^k ~ < Hq. 



3 Energy Density Perturbations 

In this section, we investigate whether the issue of structure formation is mod- 
ified by the dynamics of the BD scalar field within the BI universe model. 
This can be achieved by considering the evolution of energy density pertur- 
bations. Basically, we perturb the field equations and take the temporal 
components of these equations. The behavior of the relevant variables for 
the accelerated expansion of the universe is investigated. For this purpose, 
we introduce the notations Sg^ = h^, 5R 00 , 5u^, and ST 00 to denote the 
perturbations in g^ u , Rq , m m , and T 00 , respectively. For such perturbations, 
we use a synchronous coordinate condition given by = 0, V/i. Using 



13 



perturbation in the metric tensor defined by 

9 ' fiv = 9^f "I" fyw) 
the perturbed form of the Ricci tensor can be written as 

SRp, = X%- P ~ X P m „, (36) 

where 

The non-zero temporal component of the Ricci tensor is 

(m + 2)B S B 2 
Roo = ~ B ' + m(l - m)—. (37) 

Multiplying eq.f lT7|) by fe|fy and after some manipulation, it follows that 

2 j> p (m+l)(l-3 7B )+2(m + 2)(2o; + 3)7B 
^ 00 " + + (m + 3)(2a;o + 3) 

2 ] 2(1 + m) V ^ (m + l)V 
(m + 3) J (2wo + 3)(m + 3)> 2</> J (m + 3)0 ' 1 J 

Using the definition of the perturbed Ricci tensor given in eq.( )36|) . the per- 
turbed temporal component of the Ricci tensor is determined by 

1 B B 2 B 

6R ™ = ^^fin-4m-h 11 + 2m((l + 2m)—--)h 11 } 

+ ^[h kk -^h kk + 2(3^-^)h kk \, k = 0,2, 3. (39) 

The trace of the energy-momentum tensor is perturbed as 

T = T + 5T = T + 5p- 35P. (40) 
The action of the perturbed D' Alembertain operator on the BD field is 

Id 2 Id 2 Id 2 



B 2m d 2 x ^ B 2 d 2 y y TJ B 2 d 2 z 
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We define he following parameters for relevant perturbations 



h n = B 2m h, h kk = B l h, 5<p = A0, A«l, Sp = Ap, A < 1. (42) 

Here, h(t), A(t), and A(t) denote the perturbed gravitational, scalar and 
matter density fields, respectively. Using these parameters in independent 
field equations, we obtain 

h = A + 2(1 + ^ )4 + (A - A)g t + ^ - 3 ^ 
<p (f> (to + 3)(2w + 3) 

2(to + 2)(2o;o + 3)7b 2 _ (to + 1) _ (2-nQ 

(m + 3)(2w + 3) (m + 3) J (to + 3) [ (2w + 3) J cl 1 

-l)A0 m_1 , (43) 

1 d 2 A 9 2 A Ap(l-3 7fl ) (2-m)V c 

B 2 { d 2 y d 2 z } 0(2w o + 3) (2w + 3)' 1 J 

A - A + ^ = 0. (45) 

Here, w l is the comoving velocity of the fluid and Su\ denotes the respective 
perturbed quantity. We also introduce the notation m = For the 

discussion of structure formation (as the structure is formed in the absence of 
pressure, which can prevent gravitational collapse), we consider a pressureless 
case, i.e., 7# = 0. Thus, the parameters a and (3 turn out to be 

2 -2m , s 

P = , " = 7 r-r— r. 46 

H 1-ni' (l-m)(m + 2) v ; 

By using an infinitesimal gauge transformation, we can take the perturbed 
four- velocity as null, for which eq.( l4"5]) leads to 

A = h, (47) 

where the constant of integration is taken to be zero. We also consider a 
plane wave-like behavior of the considered perturbations, that is, 

A(x, t) = X(t) exp(-tk.x), (48) 
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where k is the wave number. Substituting eqs.( ITB|) and fiSl) into eqs. 
and (HU), it follows that 

v 4(l + q; )A (A - A) (m + 1) + 2(2u + 3) 

(l-ni)t t 2 P4 (m + 3)(2w + 3) J 

(™+V-#^]K-l)A« (49) 



(m + 3) 2w + 3 JV ' t 2 



A(4-2ni) 2A 2A AA; 2 t , 4^1 

A = — 1 1 -( — ) (l-"l)(">+2) 

t (l-ni) t 2 {i-n x ) t(i-m) B 2m H } 

\{kl + kj) t Aft MQj y^ ,(^-1) 

+ i? 2 V + (2 Wo + 3)t 2 (2o; + 3)t 2 KcAni0 ° ' 

(50) 

where p c = ^^£q ■ The expressions for p c ( m ) and in terms of ni are 
0o m 2 — 3m + 2 4n 2 2n 1 (m + 1) 



Pc(m) 



t 2 Q /{1 ' ni) 2 (1 -ni) 2 (m + 2) 2 (1 - m) 2 (m + 2) 

2 , (m + 3)m 4(1 + 6Uq) 

+ 71 — rrrd' ( 51 ) 



1 - n x (1 — ni)(m + 2) (l-ni) 5 

tf- ni) 8n 2 (l + 2m) 4n t 4a; 

c " t 2 Q [ (m + 2) 2 (l -m) 2 (1-nx) 2 (l-ni) 21 ' 1 ' 

Combining eqs. (I49p and f )50|) and neglecting the terms involving higher 
powers of t~ 2 , we have 

A + ^(A-A) + §(A-A) + ^ + ^ = 0, (53) 
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where 



1 — Tlx 



2(m + 2) + 2(2w + 3) m 2 - 3m + 2 4n 2 

O 2 (W J - [ 7— . Q w , ; . Q x J[" 



2 



(m + 3)(2o;o + 3) JL 2 (1 - m) 2 (m + 2) 
2m(m+l) 2 (m + 3)na 4(1 + w ). 

+ (1 - ni ) 2 (m + 2) ~ 1-m + (1 -ni)(m + 2) + (1 -ni) 2 ' 
(4 -2^) 4(1 + a; ) 



C 3 (wo) = Ci(a;o) 



(l-m) (l-m) 



ClW = a. w + ^- + ( '" + ! twfVif ^ + 'to 1 

1 — ni (m + 3)(2w + 3) 2u + 3 

m + 3 2u + 3 

For the solution of eq.( l53l . we take the following assumptions: 

A-A = /(t), /(t)=^ 5 , A = x t 9 , (54) 

where £ and x are arbitrary constants. The substitution of the above as- 
sumption in eq.f )53|) leads to 5 = 9 and a quadratic equation for 9 given 
by 

X 9 2 + 9[ X (C 3 - 1) + - C 3 )} + C 4 ( X - + = 0. (55) 
The roots of this equation are 

0± = ^[(x(i-C3) + £(C3-tfi)) 



± 7(X(1 - C 3 ) + £(C 3 - CO) 2 - 4 X ((C 2 - C 4 )£ - C 4 0]- (56) 

In previous papers 13 > 37 ) ; the growing modes for density perturbations are 
calculated in the asymptotic limit of |ojo| ^> 1- In a similar way, the solution 
9 in our case turns out to be 

(l-ni) x 



± l /(l-^ 2 + -^-[(m+l)(l-n 1 )(l-i)-2i]]. (57) 
X (m + 3) XX 
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This shows that in the limit \uj \ ^> 1, the asymptotic value of 6 depends on 
the power of the self-interacting potential rii and the parameter m for scale 
factors. As f3 always remains negative, n\ = will always remain positive 
and greater than 1. Also, m is a positive constant with m^l implying that 
6 + represents the growing mode for £ < x i n both cases < m < 1 and m > 1. 
Note that for m = 1, the modes we obtained match with those obtained by 
Sen and Seshadri 37 ^ and n\ — 2 corresponds to the modes obtained by 
Bertolami and Martins 13 K Thus, density perturbations can grow with time 
using the admissible negative range of uo (found in the previous section) and 
appropriate choices for the constants £, x, and n\ as well as for m in the 
asymptotic limits. It can be confirmed that the growing modes can also be 
allowed for positive BD parameter values if we take (5 > or set the arbitrary 
constants £ > X- However, the first possibility remains incompatible with 
the accelerated expansion of the model (as it can be obtained by setting 
a < 3/(m + 2) with q > 0). 

4 Summary and Discussion 

In this study, we have investigated the nature of the self-interacting potential 
and dark energy, i.e., the energy density of the BD scalar field for the BI 
universe model. For this purpose, we have considered the power law ansatz 
for the scale factors as a well as the scalar field. By taking the matter 
contents as perfect fluid, we find constraints for all the relevant parameters 
and their best-fit values. The feasible regions for these parameters are shown 
graphically for the two cases m > 1 and < m < 1. We have also discussed 
the variation in the gravitational constant and the evolution of energy density 
perturbations for the obtained model. The results are summarized as follows. 

• The obtained model is entirely based on the parameters a, (3, Uq, and 
m. Using the deceleration parameter, a > ( m + 2 ^ for an accelerated 
expanding universe model. On the basis of a, the bounds for f3 are 
found which indicate that (5 always remains negative for a > and 
m ^ 1, m > 0. By imposing positivity condition of energy density 
for both matter and scalar field, we find bounds for the BD parameter 
depending on these parameters. It is seen that for admissible ranges of 
these parameters, we have — 2 < lu < —1.4 and —1.45 < u < —1.22 
for m > 1 and < m < 1, respectively. 
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• The bounds for f3 provide a positive power law potential, which can be 
used to derive the accelerated cosmic expansion. 

• Using the best-fit values for the flatness of the universe Q m and Q^, 
we find the ranges for the EoS parameter 7^. For < m < 1, this 
parameter is constrained as —8.45 < < —3.225, and for m > 1, as 
—1.8 < 7^ < —0.23. These ranges are compatible with the ranges 
suggested by the observational data 48 ^. 

• With the help of density perturbation analysis, it is shown that 9 + 
represents the growing mode for suitable choices of parameters in both 
cases, < m < 1 and m > 1. Thus, it can be concluded that the 
modifications due to the dynamics of BD field do not disturb structure 
formation for the BI universe model. 

Note that the BD constant coupling parameter has a negative range ob- 
tained by imposing the positive energy condition and is obviously inconsistent 
with solar system experiment bounds. It can be confirmed that, if we take 
negative range of m, e.g., —2 < m < 0, then we can have positive values of 
the BD parameter, which leads to the failure of the positive energy condition 
13 ) . According to recent observations, the failure of the positive energy condi- 
tion is found to be more serious as the contribution of DE must be dominant 
as compared to matter density. This problem has widely been observed in 
the context of scalar tensor theories, e.g., 14 > 37 > 49 > 50 ). Recently, it has been 
found 51 ) that the accelerated expansion is incompatible with large values of 
the BD parameter by showing that scalar-vacuum BD equations in 5D are 
equivalent to the BD theory in 4D with a self-interacting potential and an 
effective matter field. Moreover, it has been pointed out 52 - ) that structure 
formation also remains incompatible with solar system experiment bounds. 

The solar system experiments constrain the BD parameter to be greater 
than 40,000 for a massless scalar field, but it is not true when the scalar field 
is massive. In this case, we define m ^> r _1 , where r is the experiment scale 
and m is the mass of potential as the spatial dynamics of scalar field can 
freeze for this limit. The mass scale can be taken as m^u = 10~ 27 , where 
the relevant scale is taken to be the Astronomical unit (1AU = I0 s km). 
By comparing with this mass scale, we say that the scalar field is massive 
(m ^> rriAu) or has negligible mass (m <C vtlav)- in our case, the self- 

13-2 

interacting potential is V = V c (j) p , with V c given by eq. ff25l) . t is the age 
of the universe and 2=2 > 0. If the mass of the potential term, i.e., Vss, is 
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greater than the mass scale, then the scalar field is massive, consequently, it 
could be expected that all the values of the BD parameter are admissible as 
shown in ref. 53, where the potential term is taken to be V — A0 n ; n > 0, 
which is quite similar to our case when n = ^j^- and A = V c . However, an 
exact estimation could be found using the solar system constraints of the 
Post-Newtonian parameter as shown in the literature 5i \ In that case, the 
potential is taken to be U((f>) = ^m 2 ((f) — 4>o) 2 , which is a subcase of the 
potential derived in the present study for /3 = —2 and it is shown that, for 
m > 200rriAu, all values of uj > —3/2 are admissible at a 2a confidence level. 

In our recent work 31 \ we have studied the "early- and late-time" be- 
haviors of a BI universe model using perfect, anisotropic and magnetized 
anisotropic matter contents. In the current study, we impose the condition 
a > ^2 (which is necessary for a negative deceleration parameter) and to 
determine will be the nature of the self-interacting potential V, and the en- 
ergy densities and pressures due to matter and the scalar field. By exploring 
the admissible range of the BD coupling parameter and consequently, the EoS 
parameters •jb and 7^, it is concluded that these values are well-consistent 
with the observed values. 
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